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Abstract 

A dense neutrino medium such as that inside a core-collapse supernova can experience collective flavor conversion or oscil¬ 
lations because of the neutral-current weak interaction among the neutrinos. This phenomenon has been studied in a restricted, 
stationary supernova model which possesses the (spatial) spherical symmetry about the center of the supernova and the (direc¬ 
tional) axial symmetry around the radial direction. Recently it has been shown that these spatial and directional symmetries can 
be broken spontaneously by collective neutrino oscillations. In this letter we analyze the neutrino flavor instabilities in a time- 
dependent supernova model. Our results show that collective neutrino oscillations start at approximately the same radius in both the 
stationary and time-dependent supernova models unless there exist very rapid variations in local physical conditions on timescales 
of a few microseconds or shorter. Our results also suggest that collective neutrino oscillations can vary rapidly with time in the 
regimes where they do occur which need to be studied in time-dependent supernova models. 
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1. Introduction 


Neutrinos are essential to the thermal, chemical and dynam¬ 
ical evolution of the early universe and some of the compact ob¬ 
jects such as the proto-neutron star inside a core-collapse super¬ 
nova (SN). Whenever there is a difference between the energy 
spectra and/or fluxes of the electron-flavor neutrino/antineutrino 
and other neutrino species, the flavor conversion or oscillations 
among different neutrino flavors can also have important im¬ 
pacts on nucleosynthesis and other physics inside these hot and 
dense astrophysical environments. 

In the absence of collision the flavor evolution of the neu¬ 
trino obeys the Liouville equation |[I}[3 

d,p + V ■ Vp — “i[Hvac Hyy, p], (1) 


where v is the velocity of the neutrino, p(f, x, p) is the (Wigner- 
transformed) flavor density matrices of the neutrino which de¬ 
pends on time f, position x and neutrino momentum p, Hyac 
is the standard vacuum Hamiltonian, and Hmat and Hyy are the 
matter and neutrino potentials, respectively. The neutrino po¬ 
tential in Eq. 0 takes the following form mil: 


Hyy = V2G, 


/ 


dV' 

(2nf 


(1 -V- v)[p(f,x,p')-p(f,x,p')], (2) 


where Gf is the Fermi coupling constant, and p is the density 
matrix of the antineutrino. Because the neutrino potential cou¬ 
ples neutrinos of different momenta, a dense neutrino medium 
can oscillate in a collective manner (see, e.g., CHia; see also 
Cil for a review). 

Eq. 0 poses a challenging 7-dimensional problem (not tak¬ 
ing into account the dimensions in neutrino flavors), and it has 
never been solved in its full form. In previous studies various 


simplihcations have been made so that a numerical or analytic 
solution to this equation can be found. For neutrino oscillations 
in SNe a commonly used model is the (neutrino) Bulb model 
US). In this model a spatial spherical symmetry around the 
center of the SN is imposed so that it has only one spatial di¬ 
mension. An additional directional axial symmetry around the 
radial direction is imposed to make the model self-consistent 
which reduces the number of momentum dimensions to two. 
One also imposes the time translation symmetry because the 
timescale of neutrino oscillations is much shorter than those in 
neutrino emission or dynamic evolution in SNe. However, it 
has been shown in a series of recent studies that both the spa¬ 
tial and directional symmetries can be broken spontaneously by 
collective neutrino oscillations if they are not imposed Il28l43^ 
(see also El for a short review). In both cases small devia¬ 
tions from the initial symmetric conditions are amplihed by the 
symmetry-breaking oscillation modes which can occur closer 
to the neutrino sphere than the symmetry-preserving modes do. 

It is natural to wonder if collective neutrino oscillations can 
also break the time-translation symmetry spontaneously in SNe 
Il28l . If they do, then a time-independent SN model may not ac¬ 
curately describe the neutrino oscillation phenomenon in SNe 
even though the typical timescale of the variation in the neu¬ 
trino emission is much longer than that of neutrino oscillations. 
In this letter we analyze the neutrino flavor stability in a time- 
dependent SN model which should provide some interesting in¬ 
sights to this question. 

2. Time-dependent nentrino Bulb model 

We will focus on the potential differences between the re¬ 
sults obtained from the time-dependent and stationary SN mod¬ 
els. Therefore, we will employ the time-dependent Bulb model 
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and 



Figure 1: The geometric picture of the time-dependent (neutrino) Bulb model 
for supernova. Two neutrinos emitted from the neutrino sphere of radius R with 
emission angles and and at time t - Ijc and t - I' Ic meet each other at 
radius r and time t, where / and I' are the distances by which the two neutrinos 
have traveled from the neutrino sphere to their meeting point, respectively. 


which has the same spatial spherical symmetry and the direc¬ 
tional axial symmetry as in the conventional Bulb model. Un¬ 
like the conventional stationary Bulb model, however, we will 
not assume that the emission and flavor evolution of the neutri¬ 
nos are time-independent (see Fig. [T]|. For simplicity, we will 
consider the mixing between two active flavors, the e and x fla¬ 
vors, with the latter being the linear superposition of the fi and r 
flavors. We also assume a small vacuum mixing angle 6y ^ 1. 

It is convenient to use the vacuum oscillation frequency 


\Airr\ 
' IE 


(3) 


to label the neutrino and antineutrino with energy E, where Am^ 
is the neutrino mass-squared difference, and the plus and minus 
signs apply to the neutrino and the antineutrino, respectively. 
We define reduced neutrino density matrix 


p(f; r; u, u) oc 


if m > 0, 
if m < 0 


with normalization 


trp = 1, 


(4) 


(5) 


where u - sin^ with being the emission angle of the neu¬ 
trino on the neutrino sphere (see Fig. [T]), and r is the radial dis¬ 
tance from the center of the SN. 

The equation of motion (EoM) for the (reduced) density 
matrix g can be written as 


i{df -t- Viidr)g — [Hyac + + Hyy, p], 


where 


Vu{r) 


-fE 


(7) 


is the radial velocity of the neutrino. In the weak interaction 
basis the standard vacuum Hamiltonian and the matter potential 
are 


T]U _ 
^vac ~ « ^3 “ 


rjo) 

~2 


1 0 
0 -I 


( 8 ) 


Hmat — d 


1 0 


0 

0 0 


0 0 


(9) 


respectively, where rj - +l and -1 for the normal (neutrino 
mass) hierarchy (NH, i.e. with Am^ > 0) and the inverted hi¬ 
erarchy (IH, Am^ < 0), respectively, and is the net electron 
number density. 

In this letter we assume that the number flux F'v^/v^(£’,)?«) 
of the neutrino/antineutrino in flavor a {a - e, x) is time inde¬ 
pendent ll38l . We define the distribution function of the neutrino 
emission to be 


gioj, u) > 


dE 


do) 


(Fy,+FyJ if«>0, 
-(Fy^ + FyJ if « < 0 


with normalization conditions 


^oo ^1 

du3 

Jo Jo 


y g(w,M) = 1, 

du Vf 


( 10 ) 

(lla) 

(lib) 


where 




du 

0 Jo 


du 

-Y (Fy. + FyJ, 

0 ^ 




J ^oo ^1 

dE :^(y.+y,) 
0 Jo 


(12a) 

(12b) 


du 

to Jo 2 

are the total number luminosities of the neutrino and antineu¬ 
trino (i.e. the number of neutrinos or antineutrinos emitted by 
the whole neutrino sphere per unit time), respectively. The 
opposite signs of g(oj, u) for the neutrino and antineutrino in 
Eq. ^ take into account their different contributions to the 
neutrino potential in Eq. Q. In the Bulb model the neutrino 
potential can be written as 

V2GFVf 


Hvv(f; r; u) - 


■ fd.' r ^ 

%J~oo *Jo 


47rr2 JO 

X g{cj', u’) g{t; r; w', u). 


(1 - VuVu') 


(13) 


Because collective neutrino oscillations usually occur in the 
regime R/r « 1 in the Bulb model, we will take the large- 
radius approximation ||39]| 


Vu(r) 




(6) In this approximation. 


/ (u + u' , , , 

dF , 


(14) 


(15) 


where all the primed quantities are functions of u' and a>', e.g., 
g' - g{t\ o)' ,u'\ r). 




^/2GpN^°^ (R 


AnR'^ 


(;) 


(16) 


is the strength of the neutrino potential at radius r, and 


/ r^oo 

dF' = dm' I du'. 
*J — 00 aJo 


2 


(17) 

















3. Linear regime 


In the regime where no significant flavor transformation has 
occurred, the linear flavor-stability analysis is applicable 1401 . 
In this regime the neutrino density matrices take the form 


g(t; r; ai, u) 


Qee Qxx 

1 

0 

Qee Qxx 

1 e 

2 

0 

1 

' 2 

e* -1 


(18) 


where Qee(<JJ, u) and gxx (^,«) are the probabilities for the neu¬ 
trino (or antineutrino) to be in the e and x flavors, respectively, 
and |e(f; r; at, m)| 1. Here in the spirit of flavor-stability anal¬ 

ysis we have assumed that g are approximately constant. At the 
onset of collective neutrino oscillations e grow exponentially. 
If e has strong time dependence, the time translation symmetry 
is broken spontaneously by collective neutrino oscillations. 

Keeping only the terms up to 0(e) in Eq. (|^ we obtain 


i(d, + v^dr) e (-rjio + A + C)e 

-^J(u + u')(g:^-g:jg'e'dr', (19) 

where 

c(u,m) ^ 2 “')(eee “ exx)g' dr'. (20) 

Defining 


f 


ec 7 (r; u,u) - I e(t; r; a>, u) dt 


( 21 ) 


; can rewrite Eq. in frequency space ; 


idre-sj ~ < -rjoj + (A - m) 




+ Cf e. 


J(u + u)(g^^ - g'.Jg' e^ dE'. (22) 


We note that frequency m in the above equations represents 
the temporal variation of the neutrino flavor quantum state at 
given radius r, and it shall not be confused with the vacuum 
oscillation frequency a> which is determined by the energy of 
the neutrino. 

We also note that the starting point of collective oscilla¬ 
tions is determined by the comparison of the dispersion in each 
term of the neutrino propagation Hamiltonian in Eq. Q with 
the overall strength of the neutrino potential. The spreads in 
the vacuum Hamiltonian Hyac and the neutrino potential Hw are 
dominated by variations in vacuum oscillation frequency oj and 
trajectory parameter u, respectively. Therefore, we have taken 
Vu =» 1 for these terms in Eq. ( |22l ) as in Ref. ll^ . Eor the mat¬ 
ter potential Hmat, however, the lowest order term of v^'Hmat 
in the large-radius expansion is the same for all neutrinos and 
does not suppress collective oscillations BTl . In Eq. ( |22l l we 
have included its next-order term which can suppress collective 
oscillations if a very large matter potential is present ll39l . 

Eq. ( [22| i is the same as that for neutrino oscillations in the 
stationary Bulb model except with replacement A —» A-m. The 


flavor-stability analysis of this model has been carried out in de¬ 
tails in Ref. Il40l which we shall not repeat here. The essence of 
this analysis is to find out all the collective oscillation solutions 
to Eq. ( |2^ which are of the form 

= (23) 

where QctCWj “) is independent of r, and Q.^(A,fi) is the collec¬ 
tive oscillation frequency. If 

= Im(flj^) (24) 

is positive, there exists a flavor instability, and e^ will grow ex¬ 
ponentially in terms of r which can lead to significant flavor 
transformation. If there exist multiple unstable modes, the un¬ 
stable mode with the largest exponential growth rate will 

eventually dominate. 


4. Results and Discussion 

We analyzed the flavor instabilities in the time-dependent 
Bulb model with two sets of neutrino spectra. In the first case 
we assume the same single-energy spectrum as in Ref. Il3^ 
in which all neutrinos and antineutrinos have the same vac¬ 
uum oscillation frequency wq = 0.68 km ' and the number 
fluxes are Ny^ - 1.25 x 10^^s“\ Ny^ - 8.32 x 10^^ s“' and 
Nyjy^ - 5.20 X10^^ s"'. In the second case we assume the same 
Eermi-Dirac spectra as in Refs. Il^l40ll which have degeneracy 
parameters rjy^ = 3.9, rjy^ = 2.3 and rjyjy^ - 2.1, average ener¬ 
gies = 9.4MeV, (EyJ = 13.0MeV,"<£v,/v,> = 15.8MeV, 
and luminosities = 4.1x10^'erg s“', Ly^ = 4.3x10^'erg s“*, 
Lyjy^ = 7.9 X 10^* ergs“*. In both cases we assume a neu¬ 
trino sphere of radius R = 10 km and mass-squared difference 
Am^ - -2.4 X 10“^ eV^, i.e. with an inverted neutrino mass 
hierarchy. 

In Eig. I^we plot for a few frequency modes as func¬ 
tions of neutrino potential strength g [see Eq. (dH] assuming 
that the matter density is not large enough to suppress collec¬ 
tive oscillations (i.e. v“'/l T is valid). In both cases both the 
instability region and x™” are about the same for the frequency 
modes with |td-| < 100 km *. This is not a coincidence. Com¬ 
pared to the stationary model, the time-dependent model has a 
new term in Eq. 

m (R\^ um 

— UTH- - —. (25) 

Vu \r I 2 

The first term in the above equation changes only the real part 
of the collective oscillation frequency and has no impact 
on the flavor stability. The second term depends on the neu¬ 
trino trajectory and has a spread Am ~ (R/r)'^\m\. It becomes 
important only when 

Am > Wo, (26) 

where wo ~ 1 km“* is the typical vacuum oscillation frequency 
(and also the spread of w) of supernova neutrinos with the atmo¬ 
spheric mass-squared difference. In both cases collective neu¬ 
trino oscillations occur at r ~ lOR which implies that the stabil¬ 
ity condition of the frequency modes with [nrl < 100 km“* are 
about the same. 
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Figure 2: The maximum exponential growth rates of a few Fourier modes with various frequencies m (as labeled) as functions of the strength of the neutrino 
potential ji in the neutrino Bulb model. The left panel uses a single-energy spectrum o, and the right panel a continuous spectrum I22II40I . The peaks of 
correspond to distances 107 km (left) and 149 km (right) from the center of the SN. The matter density is assumed to he not large enough to suppress collective 
neutrino oscillations. 


The above arguments can be generalized to the scenarios 
where collective oscillations occur close to the neutrino sphere 
(but not too close so that r - R ^ R) because of, e.g., spatial 
inhomogeneous oscillation modes ll3^ [33] [35] |36l or different 
angular distributions for the neutrino fluxes in different flavors 
li29ll^ . In these scenarios the spread in 'cr is of the same or¬ 
der as -uj itself, and the frequency modes with |n7| < 1 km“* a; 
(3jus) ' should have the same stability condition. These argu¬ 
ments also apply in the presence of a large matter density be¬ 
cause the comparison between Arir and coq is not affected by the 
presence of the matter potential. 

We note that there exists a causality constraint in the time- 
dependent Bulb model. Suppose that there is a temporary change 
in the neutrino fluxes on one side of the neutrino sphere which 
lasts for a time interval Af. Because it takes at least At' ~ R for 
this change to propagate throughout the proto-neutron star, the 
assumption of the spherical symmetry implies that the inequal¬ 
ity At > R must hold. Therefore, only the oscillation modes of 
frequencies 

~ (lOkmr* (30A(sr' (27) 

are allowed in the spherical Bulb model. From the above dis¬ 
cussion we conclude that there should be no significant differ¬ 
ence between the flavor stability conditions in the time-dependent 
and stationary Bulb models. For more general time-dependent 
SN models, collective neutrino oscillations should occur at ap¬ 
proximately the same radius as in the corresponding stationary 
models unless there exist very rapid variations in local physical 
conditions on the timescales of a few microseconds or shorter. 

Meanwhile, the fact that the frequency modes with |?ir| < 

1 - 100 Wo all have similar instability regions also implies that 
the time-translation symmetry can indeed be broken sponta¬ 
neously by collective neutrino oscillations in the Bulb model, 
and that neutrino oscillations can have a strong time depen¬ 
dence once collective oscillations begin. As a result, there may 
exist qualitative differences between neutrino oscillations in time- 
dependent and stationary supernova models. 
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